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Abstract. We study the zero and finite temperature Casimir force acting on a 
perfectly conducting piston with arbitrary cross section moving inside a closed 
cylinder with infinitely permeable walls. We show that at any temperature, the 
Casimir force always tends to move the piston away from the walls and towards 
its equilibrium position. In the case of rectangular piston, exact expressions for 
the Casimir force are derived. In the high temperature regime, we show that 
the leading term of the Casimir force is linear in temperature and therefore the 
Casimir force has a classical limit. Due to duality, all these result also hold for 
an infinitely permeable piston moving inside a closed cylinder with perfectly 
conducting walls. 

PACS numbers: 03.70.+k, ll.lO.Wx, 12.20.-m 



1. Introduction 

It is well known that the vacuum fluctuations of electromagnetic fields in the 
presence of boundaries give rise to Casimir force, which was shown to be attractive 
when the boundary consists of a pair of perfectly conducting parallel plates [1]. 
Since the seminal work of Casimir, many studies have been done on Casimir effect 
where different quantum fields and geometric configurations have been considered. 
Nevertheless, after 60 years of its discovery, Casimir effect is still under active re- 
search. By definition, the calculation of Casimir energy involves an infinite sum 
that needs to be regularized. However, there is still no consensus on the regular- 
ization procedures, which often lead to inconsistent results. Even for geometric 
configuration as simple as a rectangular cavity, despite zeta regularization tech- 
nique or dimensional regularization method can give finite results for the Casimir 
force acting on a wall [H [3] , some authors considered these regularization methods 
which renormalize all surface divergence terms to zero as being not physical [1]. 
Nonetheless, Geyer, Klimchitskaya and Mostepanenko 5' have recently developed 
a subtraction scheme to obtain a physically consistent Casimir force acting on a 
wall of a perfectly conducting rectangular cavity from the point of view of ther- 
modynamics. Another approach to this problem was considered by Fulling et al 
In 2004, Cavalcanti 7J proposed an alternative to this problem by adding a 
piston in the rectangular cavity. He showed that for a 2-dimensional rectangular 
piston, the Casimir force acting on the piston due to fiuctuations of a scalar field 
with Dirichlet boundary conditions is finite without renormalization and can be 
computed exactly. Since then, Casimir piston has attracted considerable interest 
IHliniliainillllllllllllllllllZllIHlITOllinilll]- it has been shown that at any 
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dimension and temperature, for either scalar field or electromagnetic field, if the 
boundary conditions assumed on all the walls are the same and is an ideal condi- 
tion, then the Casimir force is always an attractive force which tends to pull the 
piston to the nearer wall. As the length scale shrinks to the nano range, this at- 
tractive Casimir force would create undesirable effects on microelectromechanical 
and nanoelectromechanical devices known as stiction [22\ 123] , which would limit 
the functionality of the devices. As a result, there is an impelling need to search for 
scenarios that would lead to repulsive Casimir force. In the case of infinite parallel 
plates, Boyer [53^ has shown that the Casimir force acting on the plates is repulsive 
if one of the plates is perfectly conducting and the other is infinitely permeable. 
The thermal correction of the Casimir force for this configuration was considered 
in [25| I26j and it was proved that the Casimir force remains repulsive at any finite 
temperature. Another generalization of Boyer's work was considered in [27], where 
it was proved that if the plates are made of dielectric materials with nontrivial 
magnetic permeability, then for a large range of parameters, the Casimir force is 
repulsive. Similar results were obtained in |28j . For the piston geometry, Barton 
[TT] showed that the Casimir force acting on a piston made of weakly dielectric 
materials can become repulsive when the plate separation is sufficiently large. In 
this paper, we generalize the original setup of Boyer 24] to piston moving freely 
inside a closed cylinder, where the piston is perfectly conducting and the walls of 
the cylinder are infinitely permeable (see FIG.[T]). This setup has been suggested 
by Fulling et al [29| but they only considered the zero temperature Casimir force for 
rectangular piston and exact explicit formulas of the Casimir force were not given. 
In the present article, we allow the piston to have arbitrary cross section and it is 
shown that the Casimir force is free of surface divergence even without renormal- 
ization. We prove rigourously that at any temperature, the Casimir force is always 
repulsive tending to restore the piston to its equilibrium position, and the magni- 
tude of the Casimir force decreases as the piston moves towards the equilibrium 
position. In the high temperature regime, we show that the Casimir forces due to 
the blackbody radiation from the two regions separated by the piston cancel each 
other, and the leading order of the Casimir force is linear in temperature. This 
implies the existence of classical limit for the Casimir force. Due to the duality 
between electric field and magnetic field in (3+1) dimension, all the results in this 
paper also hold for an infinitely permeable piston moving freely in a closed perfectly 
conducting cylinder. 

In this paper, we use the units where ft = c = fcs = 1 in sections [5] and [31 In 
sections [4] and [H we use SI units. 

2. Cut-off dependent Casimir energy 

In this section, we first discuss the eigenfrequencies of the electromagnetic field 
for each region divided by the piston. We then give the expressions for the Casimir 
energy which was computed by exponential cut-off method in Appendix [X| 

2.1. Modes of the electromagnetic field. For a cylinder [0,L] x il, where the 
boundary surfaces {0} x O and dfl are infinitely permeable but the boundary surface 
{L} X J7 is perfectly conducting, the eigenfrequencies for the electromagnetic field 
can be written as follows. The boundary conditions are equivalent to the conditions 
n.E = and n x B = on the surfaces {0} x fl and [0, L] x dfl, and n x E = 
and n.B = on the surface {L} x ft. Here n is the unit normal vector to the 
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Figure 1: A rectangular piston moving freely inside a rectangular cylinder dividing the 
cylinder into Regions I and II. 



surface, E and B arc the electric field and magnetic field respectively. The fields 
can be expressed in terms of the vector potential A by E = — ^ and B = V x A. 
Imposing the transversality condition div A = 0, Maxwell's equations are equivalent 
to — A = 0. One can then check that for the TE modes (i.e. modes with 
Ei = 0), a set of independent solutions for A is given by Ai = 0, 



l„ — r>r>o £i 

tofj i 8X2 



Here A; G N = N U {0}. For j — 1,2,3,.. ., cpj (.^2,2^3) is a nonzero eigenfunction of 
the Laplace operator — on CI with eigenvalue coj^ j > and with Dirichlet 
boundary conditions, i.e. (PjIqq = 0. The eigenfrequency WxE.j.fe is given by 

IJ-J '^TE.j.fe — I I +^D,j- 

For TM modes (i.e. modes with Bi =0), a set of independent solutions for A is 
given by 

n (k + h) Xi , ■ , 

Ai =sin— ^ — ^J-J-^(x2,X3)e"^™-^-''*, 

TT (fc + i) TT (fc + i) Xi 1 dlhAx2,X-i) -, ^. 

A2 =—^ ^ cos — ^^-^^ -^^ e"^TMj,kt^ 

L L ujf^^^ dx2 

_7r(fc+i) _7r(fc+i)a;i 1 ai/;,(a;2, X3) 



■ cos — ^ 2 ~ ^ 



i 9X3 
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Here fc G N. For j = 1,2,3, . . ., ^'^(2^2, 2:3) is a nonconstant cigcnfunction of the 
Laplace operator on fl with eigenvalue ujff ^ > and with Neumann boundary 

conditions, i.e. = 0. ujTM.j,k satisfies an equation analogous to ([1]). 

In the case the piston has rectangular cross section, i.e., = [0, L2] x [0, L3], we 
have explicitly 



L2 J V -^3 , 



2.2. Cut-off dependent Casimir energy. The Casimir energy of the piston sys- 
tem, denoted by E^^°''^{a; Li) is a sum of the Casimir energies of Regions I and II, 
and the Casimir energy of the region outside the cylinder, i.e. 

(2) E^TW, Li) - Eir^^'^^^a) + £;^e'^"'(ii " «) + i^S^l- 

Since the Casimir energy of the exterior region -Ec^s does not contribute to the 
Casimir force acting on the piston, we need not compute it here. By definition, 
the finite temperature Casimir energy of the cylinder Eq^^'^^^ {^L) can be written 
in terms of the partition function Z associated with a canonical ensemble in the 
following way: 
(3) 

^cyU„de..(^) = -T log Z = -T log [] = ^ J] c + T J] log (l - e^/^ 

Here w runs through the set of eigenfrcquencies. The first summation on the right 
hand side of ^ is divergent. Therefore we define a cut-off dependent Casimir 
energy E^cas""^" l-^) by 

(4) i?c'as"'" (i) = ^ E + T E log (1 - e--l^ 

According to the previous subsection, we can decompose this into the sum over TE 
modes E'^^^^'^^[L) and the sum over TM modes E'^^'^^{L). In the appendix, 
we use zeta function and heat kernel techniques to compute the Casimir energies 
-^Cas™TE(-^) ^'^'^ -^Cas"TM(-^) "^^^ modcs and TM modes respectively. We 

find that at zero temperature T — Q, the contribution to the cut-off dependent 
Casimir energy of the cylinder [0, L] x 17 from the TE modes is given by 

(5) 

^/tt ' 27r ■ Stt 

fc=i j=i 

Here CD,j, j = 0, 1, 2, . . . are the heat kernel coefficients of the Laplace operator 
with Dirichlet boundary coefRcients on Q.. 7 is the Euler constant, K^{z) is the 
modified Bessel function and FPs^sq {h{s)) denotes the finite part of the function 
h{s) at s = So- Notice that the first five terms of ([5]) which contain the A — > 0+ 
divergent parts of the Casimir energy depend on L linearly. The last term tends to 
as L 00. 



^CarTE''^^^^) = — CD,oA-4 + -^CdA-' + ^Cd^2\-^ + ^ (2 log A + 7 + 2 - 2 log 2) CD A 
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For the finite temperature Casimir energy, we find tliat the cut-off dependent 
Casimir energy (|4|) of the cylinder [0,L] x due to TE modes is given by 

(6) 

Ec'^^liL) =—CD.oX-^ + -^cn,,X-' + ^cn,2\-^ + ^ (2iogA + 2 + 7 - 2iog2) cn, 
' TT y/n Zn on 

-^PP„_.(r(,w(„„-iI|:f^A-,fe 



T 
2 



=1 j=i 

00 OC 



^ ^iog(l + e-^V™^+-^>..). 

— 00 j—i 

The expressions for the TM modes contributions are simiiar with TE repiaced by 
TM and D replaced by N, where now cnj, j = 1,2,3,... are the heat kernel 
coefficients of the Laplace operator on fi with Neumann boundary coefficients, 
deleting the zero eigenvalue corresponding to constant functions. Analogous to ([5]), 
the first six terms of ([6]) depend on L linearly and the last term tends to as 
L 00. 

In general, the heat kernel coefficients cjj/jy.O: cd/n,i and cd/m,2 are known to 
be given by [30] : 

Am s{dn) ^ 

where A{il.) is the area of fl, s{dil) is the arc length of the boundary of and 

with ai the interior angle of each sharp corner of dfl and ^(7^) the curvature of 
each smooth section of d^l. For I > 3, cd/n,i can be expressed as integrals of 
functions that depend on the extrinsic and intrinsic curvatures of the boundary dQ 
and the boundary conditions imposed. In the special case where J7 is a rectangle 
[0, L2] X [0, L3], we have exphcitly 

L2L3 L2 + ^3 1 

and CDfN,i ~ for all / > 3. Here do = and 5^ — 1 since we exclude the zero 
mode from the Neumann spectrum. 

3. The Casimir force acting on the piston 

As shown by ([5]) or ([6]), the A ^ 0^ divergent terms of the Casimir energy depend 
linearly on L. This implies that A ^ 0^ divergent term of the Casimir energy of 
the piston system ^ does not depend on the piston position a. Therefore these di- 
vergences do not contribute to the Casimir force acting on the piston, which implies 
that the Casimir force acting on the piston is finite even without renormalization. 

3.1. Finite temperature Casimir force and its classical limit. From ([6]), we 
find that the Casimir force is given by 
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where i^Q^g(a) is the Umit of the Casimh' force when ii ^ cxd given by 



,8) Fs.(«)=Ti: E ^ 



Since this is obviously a positive decreasing function of a, (O shows that the Casimir 
force acting on the piston has positive sign when a < Li/2 and has negative sign 
when a > Li/2. In other words, the Casimir force always tends to restore the piston 
to the equilibrium position a = Li/2. Moreover, the magnitude of the Casimir force 
decreases as the piston moves towards its equilibrium position. 

Note that ([5]) is an exact expression for the Casimir force at any finite temper- 
ature. With the knowledge of the eigenvalues of Laplace operator with Dirichlet 
and Neumann boundary conditions on the surface fi, one can compute the Casimir 
force by this formula to any degree of accuracy. In particular, for a rectangular 
piston with cross section [0, L2] x [0, L3], we have 



00 00 00 \ C^iTy + 

(a;L2,L3)=27rT E E E 7 



^ exp zna 



'{21TY 

E E E 




exp |^27ray (2;r) 

Eq. ([5]) shows that when T — > cxd, the Casimir force is dominated by a term 
linear in T corresponding to those terms with I — 0. The remaining terms decay 
exponentially as T — *■ 00. Notice that to restore the constants h, c and fc^ into the 
expression for the Casimir force, we need to replace T by ksT/ (he) everywhere and 
multiply the overall expression by he. Therefore the high temperature expansion 
of the Casimir force is the same as the small-?i expansion of the Casimir force. A 
term with order will be accompanied with the term h}^^ . The leading term of 
the Casimir force ([5]) is linear in T implies that this leading term is independent of 
h, and the remaining terms go to zero if we formally let h goes to 0. As a result, we 
find that the Casimir force has a classical limit (or high temperature limit) given 

by 

(9) F;;!Sf"'"i(a;Li) =T V ]. 

For its small a-behavior, we derive in Appendix [A1 that 

This shows that at high temperature, when the plate separation a is small, the 
Casimir force is dominated by the term 

Interestingly, this term does not depend on the geometry of the cross section, but 
depends only on the area of the cross section. Moreover, it is equal to —3/4 times the 
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classical term of the Casimir force acting on a pair of perfectly conducting parallel 
plates. In fact, one can verify that if i^cas(a; Li) is the Casimir force when both the 
piston and the surrounding walls are perfectly conducting or infinitely permeable, 
then the Casimir force _F'cas(a; Li) when the piston is perfectly conducting but the 
surrounding walls are infinitely permeable is related to i^cas(a; Li) by 

fcas(a; Li) = 2/'cas(2a; 2Li) - Fc^^ia; Li). 

In particular, when the piston has rectangular cross section, we derive from the 
results in [20] the following alternative explicit formula for the classical term of the 
Casimir force: 



77'classical / 
^Cas 



^ ' ' ' I 167r a3 87r L2 24L3 L2 ^ ^ /fca 

I, fc2 = l fe3 = l 

^ ^ / fel=l (fe2,fc3)eZ2\{0} ^ 

X Ko (^ ^^^^^^^V (fc2£2)^ + (fc3j^3)^) I _ (a ^ Li - a) . 

We would like to remark that comparing and (fT5|) in Appendix |^ we can 
deduce that in the high temperature regime, the leading behavior of the Casimir 
force contribution from Region I alone is given by (see Appendix : 



„ 00 00 / , 

1=1 j=i ^ ^ 

~c,^dT' + ^ci,^T3 + \c2,uT-' + 0(r) + (i? - iV) . 

45 x/TT 



(10) 



This implies that as expected, the leading term of the Casimir force from Region I 
alone is the blackbody radiation term 

^A(n)T4. 

However, since the blackbody radiation from Region II would contribute a force of 
the same magnitude but opposite direction, the effect of the blackbody radiation 
on the piston cannot be observed. In fact, there is also a term proportional to 
of magnitude 

J(2x-1)T2 
6 

which has been canceled out in the final Casimir force acting on the piston. 



3.2. Zero temperature Casimir force. In the zero temperature limit, the Casimir 
force can be obtained from ([5]). When L\ — > 00, it is given by 

^^C°:f=°(«)=-^E(-1)' E I^i(2fc«^)-^E(-l)' E ^'^o(2A:aa;). 

A;— 1 <-^D,j ,i^N,j k—1 uJo,j,^^N,j 
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Table 1 
and T = OK 



The Casimir force -Fcas and the contribution from when L2 — = 0.3m 



a (m) 


Fcas (N) 


fL (n) 


io-« 


1.0238 X 104 


1.0238 X 10* 


10-* 


1.0238 X 10-12 


1.0238 X 10-12 


0.1 


9.0759 X 10-25 


1.0238 X 10-24 


0.3 


2.6683 X 10-27 


1.2640 X 10-26 



Using the same method as we derive PH)) . one finds that as a — > 0+, the leading 
behavior of the zero temperature Casimir force acting on the piston is given by 

(11) 

^oo,T=Of X _ 77r^ Cq^d + Co,jv 3Cfl(3) Ci^p + Ci^n TT C2,D + C2,N ^, Qn 
"fieri „4 ^a /z: „?. ao „2 + 'J\a ) 



Cas 



48 



960 a4 16V^ a3 

1920a4 ^ ^ 48 ^ ' 

In particular, for small plate separation, the leading term of the Casimir force is 

7^2 ^ 7 7r2 
1920a4 8 240a4 ^ 

which is —7/8 times the Casimir force acting on a pair of perfectly conducting 
infinite parallel plates, in consistent with the result of Boyer [21] for parallel plates. 
Notice that this leasing term only depends on the area but not the geometry of $7. 

In the case of rectangular piston, we have the following explicit formula for the 
zero temperature Casimir force in the Li — s- 00 limit: 



77r2 1,2^3 



1920 a" 



TT 1 

96^ 



720 Ll 



Ci^(3) 
167ri2 



1 



00 00 



2L|L| fc2=ifc3=i 



X 



X exp 



27rfc2fc3-^3 
TT{2kl + 



2^3 \ ^ 



E 



.f^0fe..3fe^UO} V(fc2i2)2 + (^3^3)2 



V(fc2i2) 



In particular, 



_ 77r2 L2L3 
Cas W -1920 a4 

in agreement with the general result (jlip . 



-Fr.,: (a) 



TT 1 

96a 



(^3^3)' 



2 +0(a°), 



4. Numerical results and Discussions 

We compute the magnitude of the Casimir force when the cross section of the 
piston is a square of dimension 0.3m x 0.3m, Li — s- 00 and the temperature is equal 
to OK, IK, 300K respectively. The results are compared to the contribution from 
the term 

pll _ 77r2 £2-^3 
^'^^ ~ 1920 a4 ' 
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Table 2: The Ceisiniir force Fcub and the contribution from the classical term FqI, 
and 7?^i----i.ll when L2 = Ls = 0.3m and T = IK. 



a (m) 


i^Cas (N) 


^classical {^J^^ 


T-iclassical, 11 /t\t\ 
^Cas y.^) 


io-« 

10-4 
0.1 
0.3 


1.0238 X 10^ 
1.0238 X 10-12 
8.1004 X 10-23 
5.9862 X 10-25 


8.9146 X 10-2 
8.9146 X 10-14 
8.1004 X 10-23 
5.9862 X 10-25 


8.9146 X 10-2 
8.9146 X 10-14 
8.9146 X 10-23 
3.3017 X 10-24 



Table 3: The Casimir force -Fcas and the contribution from the classical term Fq!^ 
and when La = is = 0.3m and T = 300K. 



a (m) 


i^Cas (N) 


77'classical (^\ 
^Cas 


TTiclassical, II /atn 
^Cas U^i 


io-« 

10-4 

0.1 
0.3 


1.0238 X 104 
2.6744 X 10-11 
2.4301 X 10-20 
1.7959 X 10-22 


2.6744 X 10 
2.6744 X 10-11 
2.4301 X 10-2" 
1.7959 X 10-22 


2.6744 X 10 
2.6744 X 10-11 
2.6744 X 10-2" 
9.9051 X 10-22 



the classical term ([9]) and the term 

7-,classical, || _ 3Cfl(3) ^ ^ ^ 

and are tabulated in Tables [1] [5] and [31 Notice that when a/L2 is small and 
aTkB/{fic) <C 1, the Casimir force is dominated by the term Fq^^, whereas if 
a/L2 is small but aTkB / {fic) ^ 1, the Casimir force is dominated by the term 
-^Cas • When a/L2 — 1, there is a considerable deviation of the Casimir force 
from the leading terms i^c^s ^classical, || ^ visible length, the Casimir force is 
negligible. However, when the separation between the piston and an opposite wall 
shrinks to lOnm, then the Casimir pressure on the piston is approximately equal 
to 1 atmospheric pressure at both zero and room temperature. 

For possible applications to nanotechnology, we also compute the Casimir pres- 
sure acting on the piston when the cross section of the piston has dimension lOOnm 
X lOOnm and the temperature is OK and 300K. The results are tabulated in Tables 
|4]and[5l At this length scale, the Casimir force is not very much affected when the 
temperature change from OK to 300K. We also find that the Casimir pressure is 
almost the same as the atmospheric pressure when a — lOnm. Reducing a at a rate 
r will result in the increase of the pressure at the rate . Since the Casimir force is 
pushing the piston outward, this will prevent the undesirable collapse of the piston 
to the opposite wall. 

FIG [2] shows the behavior of the Casimir force F^as as a function of a at T = OK 
and T = 300K. It is compared to F^^^ and ideas''''"''- We see that when a is in the 
range lOnm to 1/im, the Casimir force at both OK and 300K are dominated by the 
zero temperature parallel plate term F^^^. When a is in the range IfMo. to 0.3m, 
the force is then dominated by the classical term at T = 300K. From the graphs. 
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Table 4: The Casimir pressure Pcas and the contribution from Pq^^ when L2 — = 
lOOnm and T = OK. 



a (nm) 


Pcas (N/m2) 


pL (N/m^) 


1 


1.1375 X 10" 


1.1376 X 10'-* 


10 


1.1271 X lO'"^ 


1.1376 X 10^ 


50 


1.3248 X 10^ 


1.8202 X 10^ 


100 


2.4015 


1.1376 X 10 



Table 5: The Casimir pressure Pcas and the contribution from the classical term Pca 
and p^'^^^==''^''i'll when La = is = lOOnm and T = 300K. 



a (nm) 


Pca. (N/m2) 


P^lff ''^'^i (N/m2) 


7-,classical, 11 / o\ 

Pes (N/m2) 


1 


1.1375 X 10" 


2.9715 X 10^ 


2.9715 X 10^ 


10 


1.1273 X 10^ 


2.9641 X 10^ 


2.9715 X 10^ 


50 


1.4299 X 10^ 


1.7344 


2.3772 


100 


1.2061 X 10 


5.3876 X 10-2 


2.9715 X 10^1 



one would tend to conclude that the Casimir force is an increasing function of the 
temperature. However, this is not the case as shown by FIG [31 

5. Conclusion 

We have shown that for a perfectly conducting piston moving freely inside a 
cylinder with infinitely permeable walls, the Casimir force acting on the piston is 
a repulsive force which tends to push the piston to its equilibrium position. At 
zero temperature, when the separation a between the piston and one of its opposite 
walls is small, then the magnitude of the Casimir pressure is asymptotically equal 
to 

7 TT^hc 



8 240a4 ' 

which is the result obtained by Boyer [21] for a pair of infinite parallel plates, 
one being infinitely conducting and the other being infinitely permeable. It is —7/8 
times the zero temperature Casimir pressure acting on a pair of perfectly conducting 
parallel plates. However, at high temperature, the Casimir pressure is dominated 

by 

3 C-r(3) , „ 

when a is small. This is —3/4 times the Casimir pressure when both plates are 
perfectly conducting. It is interesting to notice the change from the ratio 7/8 at 
zero temperature to the ratio 3/4 at high temperature. We have also shown that 
at high temperature, the Casimir force acting on the piston is dominated by a term 
linear in T known as the classical term. Moreover, when a is small, the classical 
term is asymptotically equal to the area of the cross section A{il) multiply with 
the pressure (fT2|l . which is of order a~^. The correction term is of order a°. 
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L2 = 0.3m, Lg = 0.3m 





— Fcas'T = OK 




Gas 




. F^^^,T=300K 




F^^J^''^'^', T=300K 


\ %^ 
\ 




\ ^^^^ 













0.2 0.4 0.6 0.8 1 



= 0.3m, = 0.3m 





A F" 

Cas 

---Fcas'T = 300K 







0.2 0.4 0.6 0.8 1 



= 0.3m, = 0.3m 




0.002 0.004 0.006 0.008 0.01 

a(m) 
= 0.3m, Lg = 0.3m 

10-' = , ^ ^ ^ ^ ^ 1 




0.05 0.1 0.15 0.2 0.25 0.3 
a(m) 



Figure 2: The Casimir force Fcas at T = OK and T = 300K and the contribution from 
F^,^^ and Fair'"'''' plotted as a function of a. 
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a = 0.001 m, U = 0.01 m, L = 0.01 m 




0.2 0.4 0.6 

T(K) 



Figure 3: This graph shows that the Casimir force is not an increasing function of tem- 
perature. Here a = 0.001m, L2 = is =0.01m. 



For simplicity, in this paper we have assumed that the piston is perfectly con- 
ducting whereas the surrounding walls are infinitely permeable. We would like to 
consider in a future work the more general case where the piston and its surrounding 
walls are allowed to have different electric permittivity and magnetic permeability. 
It would then be interesting to determine the range of the parameters for which the 
Casimir force acting on the piston is repulsive. One can anticipate that this would 
have important applications in nanotechnology. 

Appendix A. Computations of the Casimir energy and asymptotic 

BEHAVIOR OF THE CASIMIR FORCE 

A.l. Computation of the Casimir energy. Define the following zeta functions 



oc oc 

-2s 
TEj,fc' 

J=l fe=Oj=l 
00 00 00 



Co,d(s) ^ UJjj^-, Ccylindor, Te(s) = X! X! 



fc— j—1 n— — oo 

and the heat kernels 

00 



00 



k=0 j=l k=0 

The zeta functions Cn,iv(s), Ccyiinder,TM(s), Ctm(s) and the heat kernels KQ^pf(t) 
and ifcyiindcr,TM(s) for the TM modes are defined analogously. It is well known 
that the heat kernel KQ^u{t) has the following asymptotic expansion 
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Therefore, 



cylindor.TEW = 1 + 2 V (- 1) ^^e" ^ KnA^) 
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Li V — A i-3 „ / M-2 \ , 



Consequently, the function r(s) ^cylinder, teI*) has at most shuple poles at s = ^y^, 
Z = 0, 1, 2, . . . with residues 

ReS^^3^|r(s)Ccylindor, Te(s)| = ^^C£),i- 

The A ^ 0+ behavior of the T = part of the cut-off dependent Casimir energy 
(jll can be determined as follows: 
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Here 7 is the Euler constant and FP^^_ 1 ^cylinder. te(s) is the finite part of the zeta 
function Ccyiindcr, te(s) at s = ^1/2- Since 

Ccyiindcr, TE 

-L [S) Jo 

a straightforward computation gives 

ipP,^.! Ccyiindcr, Te(s) ^^{^DA [2 - 7 - 2 log 2] - FP,=_i (r(s)Co,D (s)) } 

k=l j=l 

where K,j{z) is the modified Bessel function. Gathering the above results, we find 
that the contribution to the cut-off dependent zero temperature Casimir energy of 
the cylinder [0, L] x Vl from the TE modes is given by 

(13) 
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For the finite temperature Casimir energy, we can use the following formulas. 
On the one hand, we have (see e.g. [31]): 
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On the other hand, 
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From these, we find that the cut-off dependent Casimir energy due to TE modes is 
given by 

(16) 
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A. 2. The leading behavior of the classical term of the Casimir force at 
small plate separation. Here we compute the small-a asymptotic behavior of the 
classical term of the Casimir force ([9]) . We have 

Fair^-Ha;Li)=r E + 

= - ^ E E(-i)'-' f - t) 
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A.3. Verification of ((TOl). 

(17) 
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Appendix B. Alternative formulas for the Casimir force acting on a 

rectangular piston 

Here we present two exact formulas for the Casimir force when the piston has 
rectangular cross section. These formulas are useful when the temperature T is 
high and the separation a is small. The formulas can be used to study the behavior 
of the Casimir force when the combination aT is small and large respectively. They 
can be derived using the Chowla-Selberg formula as presented in [2D] . 
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B.2. aT > 1. 
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